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ABSTRACT 

We present a technique to identify the most probable dynamical formation scenario 
for an observed binary or triple system containing one or more merger products or, 
alternatively, to rule out the possibility of a dynamical origin. Our method relies on 
an analytic prescription for energy conservation during stellar encounters. With this, 
observations of the multiple star system containing the merger product (s) can be used 
to work backwards in order to constrain the initial orbital energies of any single, binary 
or triple systems that went into the encounter. The initial semi-major axes of the orbits 
provide an estimate for the collisional cross section and therefore the time-scale for 
the encounter to occur in its host cluster. 

We have applied our analytic prescription to observed binary and triple systems 
containing blue stragglers, in particular the triple system S1082 in M67 and the period 
distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 
and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and 
that encounters involving triples could be a significant contributor to BS populations 
in old open clusters. In general, our results suggest that encounters involving triples 
could make up a significant fraction of those dynamical interactions that result in 
stellar mergers, in particular encounters that produce multiple star systems containing 
one or more blue stragglers. 
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1 INTRODUCTION 

It has been known for some time that encounters, and 
even direct collisions, c an occur frequently between stars 
in dense stellar syste ms (|Hills fc Davll 19761 ; iHut fc VerbuntJ 
1 19831 : | Leonar dl 1 19891 . e.g.). In the cores of globular clus- 
ters (GCs), the time between collisions involving t wo single 
stars can be much shorter than the cluster lifetime (| Leonard! 
1989). The time between encounters involving binary stars 
can be considerably shorter still given their much larger cross 
sections for collision. In globular and, especially, open clus- 
ters (OCs) with high binary fractions, mergers are thought 
to occur fre quently during resonant i nteractions involving bi- 
naries (e.g. lLeonard fc Linnelf l992h What's more, collision 
products have a significant probability of undergoing more 
than one collision during a given single-binary or binary- 
binary interaction since the initial impact is expected to 
result in shock heating followed by adiabatic expansion, in- 
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creasing the cross se ction for a second collision to occur (e.g. 
iFregeau et ai]|2004l l. 

Several types of stars whose origins remain a mystery 
are speculated to be the products of stellar mergers. Blue 
stragglers (BSs) in particular are thought to be produced 
via the addition of fresh hydrogen to the cores of low-mass 
main-sequence (MS) stars. Recent evidence has shown that, 
whatever the dominant BS formation mechanism(s) operat- 
ing in both globular and open c lusters, it is likely to in some 
way depend on binary stars l|Knigge. Leigh fc Sills! 120091 : 
iMathieu fc Gellei1l2009l 'l. The currently favored mechanisms 
include collisio ns during single -binary and binary- binary en- 
counters (e.g. iLeonar df ll989l h mass transfer between the 
components of a binary system (e.g. IChen fc Han|[2008al lbl) 
and the coalescence of two stars in a close binary d ue to per- 
turbations from an orbiting tri ple companion (e.g. lEggletonl 
2006; Perets fc Fabrvckvll2009l h 

A handful of spectroscopic studies have revealed that 
in some GCs there exist BSs with masses exceeding twice 
that of the MS tu rn-off (e.g. IShara, Saffer fc Liviol Il997l ; 
iKnigge et~a l. 2008). Such massive BSs must have been 
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formed from the mergers of two or more low-mass MS stars 
since they are too massive to have been formed from mass 
transfer. In a few cases, this can also be argued for entire BS 
popula tions using photometry. For instance, IChen fc Hani 
|2008bl ) performed detailed binary evolution calculations 
to study dynamical stability during mass transfer from an 
evolving giant star onto a MS companion. Based on their 
results, it can arguably be inferred that most BSs in NGC 
188 are sufficiently bright that they probably could not have 
formed from mass transfer alone. If true, this suggests that 
most of these BSs must be the products of stellar meregers. 
Regardless of the dominant BS formation mechanism(s) op- 
erating in dense star clusters, dynamical interactions should 
play at least some role. For example, even if blue stragglers 
are formed as a result of binary evolution processes such 
as mass transfer, the progenitor binaries themselves should 
have been affected by at least one dynamical interaction over 
the course of their lifetime. 

Numerous scattering experiments have been per- 
formed to explore the outcomes of binar y-binary and, in 
particular, single-binary en c ounters (e.g. [ McMillan 1 19861 ; 
ISigurdsson fc Phinnevl 1 19931 ; iFregeau et all 12004) . Most of 
the earliest of these studies were performed in the point- 
particle limit, ignoring altogether the often non-negligible 
impli c ations of the s tars' finite sizes (e.g. iHut fc Bahcaiil 
Il983l : iMikkolal Il983l) . Later, more realistic simulations 
clearly demonstrated the importance of taking into ac- 
count the dissipative effects of tid al interactions and di- 
rect contact between st ars (e.g. iMcMillan et al.l 1 19871 ; 
IClearv fc Monaghanll990l ). As a result of the increased num- 
ber of free parameters for the encounters and the longer inte- 
gration times required to run the simulations to completion, 
few studies have been conducted to explore the outcomes 
of binary-binary encounters or interactions involving triple 
systems. 

In this paper, we introduce an analytic technique to 
constrain the most probable dynamical origin of an observed 
binary or triple system containing one or more merger prod- 
ucts. Provided the observed system is found within a mod- 
erately dense cluster environment with binary and/or triple 
fractions of at least a few percent, the probability is often 
high that it formed from a merger during an encounter in- 
volving one or more binary or triple stars. In Section [21 
we present an equation for energy conservation during in- 
dividual stellar encounters and outline the process for ap- 
plying our technique. Specifically, we present a step-by-step 
methodology to evaluate whether or not an assumed dynam- 
ical history could have realistically produced an observed 
system and describe how to determine the most probable dy- 
namical formation scenario. In Section [3l we apply our tech- 
nique to a few observed binary and triple systems thought 
to contain merger products, in particular a triple system 
that is thought to contain two BSs and the peculiar period- 
eccentricity distribution of the BS binary population in NGC 
188. We discuss the implications of our results in Section 2] 



2 METHOD 

In this section, we present a general prescription for conser- 
vation of energy during stellar encounters. We will limit the 
discussion to typical interactions thought to occur in glob- 



ular and old open clusters, although our technique can be 
generalized to any choice of parameter space. The types of 
encounters of interest in this paper will predominantly in- 
volve low-mass MS stars with relative veloc ities at infinity 
ranging from < 1 km/s to ~ 10 km/s (e.g. lLeonardlll989l ; 
Sigurdsso n fc Phinnevl [l993l ) . Our technique describes how 
to isolate the most probable dynamical formation history for 
an observed binary or triple containing one or more merger 
products by providing an estimate for the time required for a 
given interaction to occur in a realistic cluster environment. 

We begin by assuming that an observed system was 
formed directly from a dynamical interaction (or sequence 
of interactions). In this case, the observed parameters of 
the system provide the final distribution of energies for the 
system resulting from the interaction(s). After choosing an 
appropriate dynamical scenario (i.e. whether the objects in- 
volved in the interaction (s) are single, binary or triple stars), 
we can work backwards using energy conservation to con- 
strain the initial energies going into the encounter. This pro- 
vides an estimate for the initial orbital energies and therefore 
semi-major axes of any binaries or triples going into the in- 
teraction. This in turn gives the cross section for collision 
and hence the time required for the hypothesized interac- 
tion^) to occur. 

Since the formation event must have happened in the 
last tbs years, where tbs is the lifetime of the merger prod- 
uct, a formation scenario is likely only if the derived en- 
counter time-scale is shorter than the lifetime of the merger 
product (s). Conversely, if the derived encounter time-scale 
is longer than the lifetime of the merger product (s), then 
that dynamical formation scenario is unlikely to have oc- 
curred in the last tbs years. In general, the shorter the de- 
rived encounter time-scale, the more likely it is that one or 
more such encounters actually took place within the lifetime 
of the merger product (s). Finally, if the derived encounter 
time-scale is longer than tbs for every possible dynamical 
formation scenario, then a dynamical origin is altogether un- 
likely for an observed multiple star system containing one 
or more BSs. Either that, or the encounter time-scales must 
have been shorter in the recent past (or, equivalently, the 
central cluster density must have been higher). 



2.1 Conservation of Energy 

Consider an encounter in which at least one of the two bod- 
ies involved is a binary or triple star. Though a complex 
exchange of energies occurs, energy must ultimately be con- 
served in any dynamical interaction. The total energy that 
goes into the encounter must therefore be equal to the total 
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energy contained in the remaining configuration: 

Ni N f 

Mi (Ii , Wj ) + ]T W i ( m * . X i . Z * > r * ) 

i i 

^ Ui(m,i,X iy Zi,n) + ^ Tj(M c i, Mj)+ 

i j 
Mi N f 

^2 e k(fJ,k, Mfc, afe) = ^ (i) 
^ Wu(mii, Xu, Za, To) + y~] Uii(mu, Xa, Zu,m) 

ii ii 

+ Tjj(Mjj) + £fcfc(/Jfcfc, Mfcfc, <2fcfc ) + A, 
M fcfc 

where A?i, Si and Mi are the total number of stars, objects 
(single, binary, triple or even quadruple stars) and orbits, 
respectively, that went into the encounter. Similarly, Nf, Sf 
and Mf are the total number of stars, objects and orbits 
remaining after the encounter. 

We let fij = Qi(Ii,LUi) represent the bulk rotational 
kinetic energy in star i, which is a function of the star's 
moment of inertia Ii and angular rotation rate u^: 



Qi(Ii,UJi) = -jli^i 



(2) 



The moment of inertia is in turn a function of the density 
profile within a star, which changes along with its internal 
structure and compos ition as the star evolves. The moment 
of inertia is given by (|Claret fc Gimenezlll989l ): 



I = ^ J p(r)r 4 dr 



(3) 



= /3 2 mR 2 , 



where j3 is the radius of gyration. For example, a typical 1 
Mq star in an old open cluster having an age of ~ 6 Gyrs 
has p = 0.241. For comparison, a 1 Mq star in a typical 
Milky Way GC haying a n age of ~ 11 Gyrs has /3 = 0.357 
l|Claret fc Gimene3 ll98gj) . A large spread of rotation speeds 
have been observed for MS stars in open and globular clus- 
ters, with measured val ues ranging from ~ 1 - 20 km s _1 
l|Mathieu fc Gelled I2009T ). A 1 M star having a radius of 1 
R.0 and a rotation speed of 2 km s _1 with (3 = 0.241 has 
Q, ~ 10 36 Joules. 

We let Wi = Wi(mi,Xi, Zi,Ti) represent the gravita- 
tional binding energy of star i, where m% is the star's mass, 
Xi is its initial hydrogen mass fraction, Zi is its initial metal- 
licity and Tj is its age. For a spherical mass with a density 
distribution pi(r), the gravitational binding energy is given 
by: 



W l (m l ,X i , Zi,n) = G j p l (r)r 4 dr 



(4) 



-S(mi,Xi,Zi,Ti)G 



R, 



where the parameter S is chosen to reflect the structure of 
the star and is therefore a function of its mass, age and 
chemical composition. For instance, a typical 1 Mq star with 
(X, Z) — (0.70, 0.02) in an old open cluster with an age 
of ~ 6 Gyrs has S = 1.892. For comparison, an older but 



otherwise identical star in a typic al Milky Way GC with a n 
age of ~ 11 Gyrs has S = 6.337 (| Claret fc Gimenej|l989h . 
Roughly regardless of age, this gives \W\ ~ 10 42 Joules for 
a 1 Mq star with a radius of 1 Rq. 

We let Ui — Ui(rrii, Xi, Zi,Ti) represent the total inter- 
nal energy contained in star i (i.e. the star's thermal energy 
arising from the random motions of its particles) . By solving 
the equations of stellar structure, the total internal energy 
of a purely isolated single star is uniquely determined by 
its mass, initial composition and age. Stars are made up 
of a more or less virialized fluid so that, ignoring magnetic 
fields, the gravitational binding energy of a star in hydro- 
static equilibrium is a bout twice its internal thermal energy 
l|Chandrasekharll 19391) . Using this version of the virial theo- 
rem, a 1 Mq star having a radius of 1 R© has U ~ 5 x 10 41 
Joules. 

The total translational kinetic energy of object j (single, 
binary, triple, etc. star) is represented by Tj- 



T j (M al ,M j ) = ^M j v], 



(5) 



where Mj is the total mass of the object, and Vj is its bulk 
translational speed. The translational velocities of stars in 
clusters for which energy equipartition has been achieved as 
a result of two-body relaxation obey a Maxwell-Boltzmann 
distribution, with the heavie st stars typica lly having the low- 
est velocities and vice versa <|Spitzerll 19871 ). According to the 
virial theorem, the root-mean-square velocity v rma of the 
distribution depends on the total mass of the cluster M c i. It 
follows that the velocities of stars in a fully relaxed cluster 
are approximately determined by their mass and the total 
mass of the cluster. Assuming that the typical velocity of a 
< m > star is roughly equal to v rms , where < m > is the 
average stellar mass in the cluster, energy equipartition can 
be invoked in some clusters to approximate the translational 
kinetic energy, and hence velocity, of a star or binary of mass 



Mi 



/< m >\V 2 
V Mi ) 



(6) 



We note that Equation [6] can only be applied in clusters for 
which the half-mass relaxation time is much shorter than 
the cluster lifetime. At the same time, the tidal truncation 
of the velocity distribution must not be significant. We will 
return to this in Section \2. 21 

For a 1 Mq star with a translational speed of 1 km 
s _1 (typical of stars in old open clusters), we find T ~ 10 36 
Joules. For comparison, the same star traveling at a speed of 
10 km s -1 (typical of stars in globular clusters) has T ~ 10 38 
Joules. We can put this into perspective by equating Equa- 
tion [5] with Equation |4j which shows that a direct collision 
between two 1 Mq stars would require an impact velocity 
of ~ 1000 km s" 1 in order to completely unbind the merger 
remnant. 

The total orbital energy of orbit k is denoted , and is 
given by: 

GfikMk 



e fe (^ fe , M fe , a,k) 



2a k 



(7) 



where fit ~ mi rail (mi + 7712) is the reduced mass of the 
orbit, Mfc = mi + m.2 is the total mass and a k is the orbital 
semi-major axis. A binary composed of two 1 Mq stars with 
a period of 1000 days has lei = 10 38 Joules. For comparison, 
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an otherwise identical binary with a period of 1 day has 
|e| = 10 40 Joules. Since most stable triples are observed to 
have outer perio d s of ~ 1000 days with clo se inner binaries 
|Tokovininlll997l ; IPerets fc Fabrvcky||2009h . it follows that 
the total orbital energies of stable triples will be dominated 
by the orbital energies of their inner binaries. 

After the encounter occurs, the energies that collectively 
define the state of the newly formed system include the rota- 
tional kinetic wm)), thermal (f/«(m«, X i{ , Zu,Tu)) 
and gravitational binding (Wu{mu, Xa, Zn, Tu)) energies 
of the stars that are left-over, as well as the total trans- 
lational kinetic (Tjj(mjj)) and orbital (e kk ([i kk , M kk , a kk )) 
energies of any left-over stars, binaries or triples. The inter- 
nal and gravitational binding energies of the left-over stars 
once again depend on their mass, composition and age. The 
mass, composition and evolutionary status of a star will de- 
cide how it responds to tidal interactions (and therefore how 
much tidal energy is deposited) since these are the princi- 
pal factors that determine its internal structure, in partic- 
ular whether or not its envelope is radiative or convective 
( Podsiadlowski 1996). Finally, we let A represent the energy 
lost from the system due to radiation and mass loss. We do 
not expect this term to be significant for the majority of the 
encounters we will consider since the rate of mass-loss from 
low-mass MS stars is small and the time-scales under con- 
sideration are relatively short compared to the lifetimes of 
the stars. Moreover, the velocity dispersions characteristic 
of the clusters we will consider are sufficiently small that we 
expect most collisions to have a relatively low impact veloc- 
ity. Though significant mass loss can occur for high impact 
collisions as often occur in the Galactic centre, mass loss is 
only < 5% for the low impact velocity c ollisions expecte d 
to occur in open and globular clusters (|Sills et alj l200lt ). 
Therefore, we will henceforth assume A ~ 0. 

Depending on the initial parameters of the encounter, 
one or more terms in Equation [T] can often be neglected. For 
example, given that the rotational energies of typical MS 
stars in both open and globular clusters are several orders 
of magnitude smaller than the other terms in Equation [T] 
we can neglect these terms for the types of encounters of 
interest. Even a 1 Mq star rotating at a rate of 100 km s _1 
has a rotational energy of only ~ 10 39 Joules. This is consid- 
erably higher than even the highest rotat ion rates observed 
for b oth BSs and normal MS stars (e.g. iMathieu fc Gellerl 
2009). From this, we expect typical rotational energies to be 
significantly smaller than the orbital energies of even mod- 
erately hard binaries. Furthermore, Equation [5] can be com- 
bined with Equation [5] and equated to Equation [7] in order 
to define the hard-soft boundary for a given cluster. If a bi- 
nary is hard, its initial orbital energy will outweigh its initial 
translational kinetic energy during typical stellar encounters 
(HeEEic 197§j). We do not expect soft bin aries to surviv e for 
very long in dynamically-active clusters (|Heggidll975l ) (i.e. 
most binaries are hard), so that the translational kinetic en- 
ergies of the stars, binaries and triples typically found in 
OCs can usually be neglected when applying Equation [T] 
This suggests Tj — 0. However, the final translational kinetic 
energies of the stars, binaries and triples left-over should not 
be neglected for encounter outcomes in which one or more 
stars are ejected with very high velocities. This can leave the 
remaining stars in a much more tightly bound configuration 
since stars ejected with high velocities carry off significant 



amounts of energy. Finally, provided there are no very high 
impact collisions and tides dissipate a negligible amount of 
energy from the system, we also expect that: 

Ni Ni N f N f 

J2u,+J2 w *~J2 u " + J2 w "- ( 8 ) 

i i an 

We can simplify our energy conservation prescription 
considerably for the majority of encounters occurring in old 
OCs. Neglecting the rotational kinetic energies of the stars, 
Equation [T] becomes: 

Mi M f 

^ e k (fj, h ,M k ,ah) = 22 £fcfc(/ifcfc, M kk , a kk )+ 

Sj kk (9) 

n 

where we have assumed Tj = and A m is the amount of 
energy deposited in any merger products formed during the 
encounter. In other words, the term A m provides the re- 
quired correction to Equation [S] resulting from internal en- 
ergy being deposited in the merger product (s) as a result 
of collisions and tides. If no mergers occur or all of the or- 
bital energy of the merging binary is imparted to the other 
interacting stars, A m = 0. 

Both the total linear and angular momenta must also 
be conserved during any dynamical interaction. This pro- 
vides two additional constraints that must also be satisfied, 
however the total linear and angular momenta are vector 
quantities that depend on the angle of approach as well as 
the relative orientations of the objects. As a result, there 
are more free parameters to fit when trying to constrain the 
initial parameters of an encounter using conservation of lin- 
ear and angular momentum. Therefore, it is considerably 
more difficult to extract information pertaining to the ini- 
tial orbital parameters of an encounter using conservation 
of momentum than it is with conservation of energy. 

We have assumed that exchange interactions do not oc- 
cur in applying our technique. Clearly, this assumption could 
be invalid for some systems. If this is the case, then our 
method is also invalid. This suggests that our technique is 
ideally suited to clusters for which the encounter time-scales 
are comparable to the lifetime of a typical merger prod- 
uct. This maximizes the probability that BSs do not expe- 
rience any subsequent dynamical interactions after they are 
formed. Of course, most clusters of interest are unlikely to 
satisfy this criterion. We can assess the probability that an 
exchange encounter has occurred for a particular system by 
calculating the different encounter time-scales and compar- 
ing to the lifetime of the merger product. If the time-scales 
are sufficiently short, then the possibility of an exchange in- 
teraction having occurred after the system's formation must 
be properly addressed. We will discuss this further in Sec- 
tion n 

2.2 Generalized Approach 

In this section, we outline a step-by-step methodology to 
constrain realistic dynamical formation scenarios that could 
have resulted in the production of an observed stellar system 
containing one or more merger products. These steps are: 
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(i) We must first find qualitative constraints for the sys- 
tem's dynamicai history and, in so doing, converge on the 
most probabie formation scenario. The choice of formation 
history should be guided by the observed properties of the 
binary or triple system containing the merger product (s). 
The following guidelines can be applied to find the most 
probable scenario. 

First, the analytic rates for single-single (1+1), single- 
binary (1+2), single-triple (1+3), binary-binary (2+2), 
binary-triple (2+3) and triple-triple (3+3) encounters can 
be compared to obtain a rough guide as to which of these 
encounter types will dominate in a given cluster. The total 
rate of encounte rs of a given ty pe in a cluster core is well ap- 
proximated by (jLeonar dl ll989r ) (see Appendix [£J for a more 
generalized form for this equation): 



r A^QTi[)^7gj , (D re ; !rms )^ re / rms , 



(10) 



where iVo is the number of single, binary or triple stars in 
the core and no is the mean stellar, binary or triple num- 
ber density in the core. The gravitationally-focu sed cross 
sectio n for collision a g f is given by Equation 6 of iLeonardl 
Gravitationally-focused cross sections for the various 
encounter types are provided in Appendix [X] along with the 
values assumed for their pericenters. 

In general, the number of single, binary and triple stars 
are given by, respectively, (1 - f b - ft)N c , f b N c and f t N c , 
where iV e = 2/37TOq r;j is the total number of objects in the 
core (jLeonardl [l989) . is the fraction of objects that are 
binaries and ft is the fraction of objects that are triples. As- 
suming for simplicity that v re i is roughly equal for all types 
of encounters, the rates for two types of encounters can be 
compared to find the parameter space for which one type of 
encounter will dominate over another. These relations can be 
plotted in the f b — ft plane in order to partition the param- 
eter space for which each of the various encounter types will 
occur with the greatest frequency, as illustrated in Figure [1] 
Given a cluster's binary and triple fractions, this provides 
a simple means of finding the type of encounter that will 
occur with the greatest frequency. Our results are in rough 
agreement with those of ISieurdsson fc Phinnevl (|l993r i who 
found that single-binary interactions dominate over binary- 
binary interactions in clusters having core binary fractions 
fb < 0.1, and may dominate for fb up to 0.25-0.5 in some 
cases. 

Second, the masses of the components of an observed bi- 
nary or triple system provide a lower limit for the number of 
stars that could have gone into its formation. The minimum 
number of stars that must have merged to form a given colli- 
sion product, labelled N m i n ^, is equal to the integer nearest 
to and larger than the quantity m rem /mTO, where m rem 
is the mass of the merger remnant and mro corresponds 
to the mass of the main-sequence turn-off (MSTO). This 
assumes that mro has not changed significantly since the 
dynamical formation of the system, which should be valid 
provided the merger products are significantly more massive 
than the turn-off so that their lifetimes are relatively short. 
The number of stars that went into an encounter Ni must 
therefore satisfy: 



0.4 0.6 
Binary Fraction 



Figure 1. Plot showing the parameter space in the — ft plane 
for which each of the various encounter types dominate. Bound- 
aries between regions are indicated by solid lines, each segment 
of which is obtained by equating two particular encounter rates 
using Equation I10| and the r elevant cross sections derived using 
Equation 6 of Leonard ( 198§). We assume at = 5ab and ab = 90-R 
in obtaining the relations between encounter rates. This is a rea- 
sonable choice for the ratio between the average binary and triple 
geometric cross sections given that the ratio between the outer 
and inner orbital semi-major axes of triples must be relatively 
large (by a factor of > 10) in order for them to be stable l|Mardlind 
120011) . 



where Mo is the number of normal stars (i.e. not formed from 
mergers) and M denotes the number of merger products. 

Third, an estimate for the average stellar mass, and 
therefore the masses of typical stars expected to un- 
dergo encounters, should be guided by a realistic stellar 
mass function for the host cluster. First of all, observa- 
tions have shown that a significant depletion of low-mass 
stars occurs in dynamically evolved clusters (i.e. those for 
which t r h tage, where t r h is the half-mass relaxation 
time and t age is the cluster age) since they are prefer- 
enti ally ejected from the cluster during close encounters 



2005 



von Hippcl & Sarajcdini 1998; Bonatto, Bica & Santos 



De Marchi. Paresce fc Portegies ZwartJ 120101 ). This 



+ M , 



(11) 



suggests that OCs and the least massive Milky Way GCs 
should have stellar mass functions that, at least in their 
central cluster regions, appear eroded at the low-mass end. 
From this, we expect < m >< mro- Finally, most of the 
stars involved in dynamical encounters should have masses 
close to or even slightly greater than the average stellar mass 
< m >. This is because gravitational focusing is strongest 
for the most massive objects, contributing to a shorter en- 
counter time-scale. 

Finally, the most likely formation scenario will, strictly 
speaking, minimize the total or cumulative encounter time- 
scale. However, the total time-scale required for a dynamical 
formation scenario that involves more than one encounter 
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will typically be dominated by the second encounter. This 
is because, after an initial encounter has occurred between 
any two suitable objects to form a new stellar configuration, 
the time-scale for a second encounter to occur is given by 
the time required for the product of the initial encounter to 
experience a subsequent encounter. This increases the en- 
counter time-scale by a factor No. In other words, Equa- 
tion [10] provides the time required for any two of the spec- 
ified objects to experience an encounter (two binaries, a bi- 
nary and a triple, etc.). It can be multiplied by iVo to ob- 
tain an estimate for the time required for a specific object 
to experience an encounter. For example, to find the time 
required for a specific binary to experience an encounter 
with another single, binary or triple star, we must multiply 
by a factor fbN a . It follows that the timescale for multiple 
encounters to occur is considerably longer than any of the 
single encounter timescales. Therefore, unless the number 
of either binary or triple stars is very low, scenarios involv- 
ing the fewest number of encounters are generally preferred 
since this tends to minimize the total time required for the 
encounter(s) to occur in a realistic cluster environment. 

The times between the various types of encounters can be 
derived using Equation [TTJ and have been provided in Ap- 
pendix [S] These time-scales can be multiplied by (1 — /j, — 
ft)N c , fbN c or ftN c to find the time required for a particu- 
lar single, binary or triple star, respectively, to encounter an- 
other object. This yields tim e-scales that are in rough ag ree- 
ment with Equation 8-125 of lBinnev fc Tremaind (|l987lT . As 
an example, we can multiply Equation IA8I by fbN c to find 
the time for a particular binary to experience an encounter 
with another binary. This gives: 

V n ) (12) 

\5km/s) \< m >/ \ a b / years ' 

Provided the derived encounter time-scale is shorter than 
the lifetime of the merger product(s), this suggests that the 
encounter scenario in question could be realistic and is there- 
fore a candidate formation history. Conversely, if the derived 
encounter time-scale is longer than the lifetime of the merger 
product(s), then the dynamical formation history is unlikely 
to have actually occurred. In general, the shorter the derived 
encounter time-scale, the more likely it is that one or more 
such encounters actually took place within the lifetime of 
the merger product (s). 

The preceding guidelines specify a narrow range of allowed 
formation scenarios. In particular, they constrain the num- 
ber of stars involved in the encounter(s), the types of objects 
(i.e. single, binary or triple stars) involved, and the number 
of encounters that took place. In most cases, these guidelines 
will converge on a single qualitative formation history that 
is unique up to the possible initial distribution of energies 
that describe suitable interactions. 

(ii) Next, we must assign an approximate value based on 
the observations to every parameter possible in Equation [TJ 
Nearly all of the required information pertaining to the fi- 
nal distribution of energies in Equation [JJ can be found from 
spectroscopy alone, though repeated measurements spread 
out over a sufficiently long timeline will typically be required 
to obtain orbital solutions and to detect outer triple compan- 
ions whenever they are present. This gives the final orbital 



energies euk as well as the gravitational binding energies Wu 
of the stars according to Equation[7]and Equation[4] respec- 
tively. Since merger products have be en shown to typ ically 
be in hydrodynamic equilibrium (e.g. ISills et al-lfeOQlf ). the 
stars' internal energies Uu can then be approximated using 
the virial theorem. Alternatively, stellar models can be used 
together with photometry. 

The measured broadening of spectral lines gives an es- 
timate for the stars' rotation speeds (although there is a 
strong dependence on the angle of inclination of the stars' 
axis of rotation relative to the line of sight), which in turn 
provides their rotational kinetic energies flu according to 
Equation [2] In conjunction with proper motions, radial ve- 
locity measurements also provide an estimate for the sys- 
temic velocity of the final stellar configuration relative to 
the cluster mean, which in turn gives its translational kinetic 
energy Tjj according to Equation [5] Finally, the cluster ve- 
locity dispersion provides an estimate of the relative velocity 
at infinity for a typical encounter, which in turn decides the 
initial translational kinetic energies Tj of stars or binaries 
involved in the encounter. Under the assumption of energy 
equipartition for both single and binary stars, the initial ve- 
locities of the impactors can be approximated by Equation[6] 
The assumption of energy equipartition should be valid in 
clusters for which the half-mass relaxation time is consider- 
ably shorter t han the cluster lifetime, and this is the case 
for most GCs (|Harrisill996l ; iDe Angeli et all 120051 ) . On the 
other hand, this assumption is likely invalid for most open 
clusters. This is because the tidal truncation of the veloc- 
ity distributions are significant in OCs since they are much 
less centrally concentrated. In this case, the initial veloci- 
ties of the impactors can be approximated from the velocity 
dispersion, which is nearly independent of mass. 

(iii) We can now obtain quantitative constraints for the 
initial encounter(s). Once we have decided on a qualitative 
encounter scenario that could have produced the observed 
merger product (s), we can estimate the orbital energies of 
the initial binaries or triples going into the encounter us- 
ing Equation [5] This gives us an equation that relates the 
initial orbital semi-major axes of all orbits going into the 
encounter (to each other) . If only one orbit goes into the en- 
counter, then we can solve for it explicitly. From this, we can 
constrain the initial collisional cross sections for realistic en- 
counters. If the derived cross section is significantly smaller 
than the average semi-major axes of all binaries and/or 
triples in the cluster, then we can infer that the time re- 
quired for the encounter to occur is significantly longer than 
the corresponding time-scale given in AppendixfX] If we use 
our derived cross section found from Equation [9] instead of 
this average semi-major axis, this should give us an idea of 
the time required for that particular type of encounter to 
occur. Strictly speaking, however, this is only a rough ap- 
proximation since these are typical time-scales found using 
the average period (or, equivalently, semi-major axis and 
hence cross section). 

If we find that the derived encounter time-scale is longer 
than the lifetime of a typical merger product, then the cho- 
sen formation history is unlikely to have actually occurred. 
In this case, it becomes necessary to re-evaluate the possi- 
ble dynamical formation histories of the observed binary or 
triple system, choosing the next most likely qualitative sce- 
nario for further quantitative analysis. These steps can be 
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repeated until either a suitable formation history is found or 
the list of possibilities is exhausted so that the only remain- 
ing conclusion is that the observed binary or triple system 
is unlikely to have a dynamical origin. 



3 RESULTS 

In this section, we apply our technique to two particular 
cases of observed binaries and triples containing merger 
products. The first is an observed triple system in the 
old open cluster M67 that is thought to contai n two BSs 
|van den Berg et al.l [200ll ; ISandauist et all 120031 ). The sec- 
ond is the period-eccentricity distribution of the BS binary 
population of the O C NGC 188, which bears a remarkable 
resemblance to M67 (jMathieu fc Gellerl [20091 ) . After deter- 
mining the most probable qualitative formation scenarios, 
we obtain quantitative constraints for suitable initial condi- 
tions that could have produced the observed orbital param- 
eters. 



3.1 The Case of S1082 

S1082 is believed to be a t riple system in the old OC M67 
jvan den Berg et alJlioOll ; ISandauist et al.ll2003l ). The ob- 
servations suggest that a distant triple companion orbits 
a close binary containing a BS and another peculiar star. 
The companion to the BS has a photometric appearance 
that puts it close to the MSTO in the CMD and yet, cu- 
riously, its derived mass is significantly greater than that 
of the turn-off. The outer companion is a BS in its own 
right, so that S1082 is thought to be composed of two BSs. 
Although both the inner and outer components of this sus- 
pected triple have systemic velocites that suggest they are 
both cluster members, it is important to note that there is 
no direct evidence prov ing a dynamical link between the two 
ijSandquist et al.ll2003l ). Assuming for the time being that a 
dynamical link does exist, we can apply the procedure out- 
lined in Section (231 to the case of S1082: 

(i) Before applying our technique, it is important to con- 
vince ourselves that a dynamical origin is possible for the 
observed system. This is certainly the case for S1082 since 
no known BS formation mechanism could have produced 
the observed stellar configuration without at least some help 
from dynamical interactions. The first step of our procedure 
is to find qualitative constraints and, in so doing, isolate the 
most probable encounter scenario. 

First, we need to know the cluster binary and triple frac- 
tions in order to use Figure [T] to find the encounter type 
occurring with the greatest frequency. From thi s, we find 
that 2 +2 encounters presently dominate in M67. iFan et all 
|l996h showed that observations of M67 are consistent with 
a cluster multiple star fraction ~ 50%. More recent studies 
report a l ower limit for fb + ft that is consistent with th eir re- 
sults (e.g. lLathamll2005l ; lDavenport fe Sandquistll2010l ). Ra- 
dial velocity surveys and simulations of dynamical interac- 
tions suggest that old OCs like M67 are likely to hos t a num- 
ber of triples with outer periods < 4000 days (e.g. iLathaml 
l2005l ; llvanova et al.ll2008l). We as sume f t /f b ~ 0.1 and, us- 
ing the result of Fan et all (| 19961 ). this gives ft, ~ 0.45 and 
ft ~ 0.05. Our assumed ratio ft/fb is slightly lower than 



found for the field, or ft/fb ~ 0.2 (Eggle ton fc Tokovininl 
l200Sfl . 

Second, we need to constrain the number of stars that 
went into the encounter. The mass of the MSTO in M67 is 
estimated to be ~ 1.3 M dMathieu fc Gellerll2009l). The to- 
tal mass of S1082 is ~ 5.8 M m |Sandquist et al.ll2003l ). From 
Equation 111! its formation must therefore have involved at 
least 5 stars. 

At this point, we can conclude that a single 2+2 encounter 
could not have produced S1082 since this scenario involves 
only 4 stars and we know that at least 5 stars are needed. 
We must therefore consider either a single 2+3 or 3+3 en- 
counter, or a multiple encounter scenario. In order to isolate 
the most probable of these possibilities, we must calculate 
and compare their encounter time-scales. To do this, we re- 
quire estimates from the observations for a few additional 
clus ter parameters. For M67, the core radius is r e ~ 1.23 
pc (|Bonatto. Bica fc Santos! 120051 : IG iersz. Heggie fc Hurlevl 
2008). From this and the central velocity dispersion, we can 
ca lculate the central mass d ensity using Equation 4- 124b 
of lBinnev fc Tremaind <| 19871 ). which gives po ~ 10 1 ' 9 Mq 



pc 



" 3 . The central stellar number density can then be ap- 



proximated according to: 



n = 



Po 



M 



< m> L 



(13) 



where M/L is the cluster mass-to- light ratio and should be 
aroun d 1.5 for an OC as old as M67 (e.g. iDe Griis et al.l 
2008). We take f t ~ 0.05 and assume most stable triples 
have outer periods of P ~ 1000 days so that at ~ 3 AU 
(assuming all three stars have a mass of 1 Mq). We also 
take ft ~ 0.45 and assume an average binary period of 
P ~ 100 days so that a b ~ 0.6 AU (assuming both compo- 
nents have a mass of 1 Mq). We assume an average stellar 
mass of < m >~ 1.0 M q, which is in reaso nable agreement 
with the observations l)Girard et al.lfl989T ). Assuming that 
the average mass of merger rem nants is equal to 2 < m > 
and extrapolating the results of ISills et al.l (|200ll ) for solar 
metallicity and more massive parent stars, we will assume 
that the typical lifetime of a merger product is tbs ~ 1-5 
Gyrs. 

A comparison of the relevant encounter time-scales sug- 
gests that the most probable dynamical formation scenario 
for S1082 is a single 2+3 encounter, although a single 3+3 
encounter is almost equally as probable. Given our assump- 
tions, we find 8.9 x 10 8 years and 3.3 x 10 9 years for the 
times between 2+3 and 3+3 encounters, respectively, in the 
core of M67. From this, we expect approximately two and 
zero 2+3 and 3+3 encounters, respectively, to have occurred 
within the last tbs years. 

From Equation 1121 we find that the time for a particu- 
lar binary to encounter another binary is 7.3 x 10 10 years. 
Similarly, the time for a particular quadruple to encounter 
another binary is 8.9 x 10 9 years (using Equation [12] and 
assuming the quadruple has a mass of 4 < m > and its ge- 
ometric cross section is twice as large as the average outer 
semi-major axis of triples). Since these time-scales are con- 
siderably longer than the cluster lifetime, this suggests that 
a formation scenario for S1082 involving back-to-back 2+2 
encounters is unlikely, even if the second encounter occurred 
sufficiently soon after the first that all four stars comprising 
the initial pair of interacting binaries are still gravitation- 



8 Nathan Leigh and Alison Sills 



ally bound. If we replace one of the 3 binaries involved in 
this scenario with a triple system, the total encounter time 
remains longer than the cluster lifetime. 

(ii) Before more quantitative constraints can be found, we 
must refer to the literature in order to obtain estimates for 
every term in Equation [9] The observations suggest that a 
binary system composed of a 2.52 Mq blue straggler (com- 
ponent Aa) and a 1.58 M q MS star (component Ab ) has a 
period of P ~ 1.068 days (van den Berg et al.ll200lh . There 
is also evidence for a 1.7 Mq blue straggler companion (com- 
ponent B) that ac ts as a stable outer tr iple with a period of 
P ~ 1188.5 days (jSandauist et al.ll2003h . From this, we can 
calculate the final orbital energies (e^ and e_e) of the triple: 

M + \e B \ ~ M ~ 10 41 J, (14) 

where we have used the fact that \ea\ 2> |es| since 
\cb\ ~ 10 39 J. It is important to note that these pecu- 
liar stars are often underluminous, so the inferred mass 
of the tertiary companio n should be taken with caution 
l|van den Berg et alJIioOlh . 

The systemic radial velocit y of the system is 33.76 ± 0.12 
km/s (jSandauist et al.l 120031 ). Although this only provides 
us with an estimate for the systemic velocity of S1082 in 
one dimension, it is consistent with the cluster mean ve- 
locity of 33.5 km/s. Therefore, the available knowledge for 
the systemic velocity of S1082 is consistent with its final 
translational kinetic energy being negligible. From this, we 
take Ta,b ~ 0. The central velocity dispersion in M67 is 
only 0.5 km s" 1 (|Mathieu fc Latham! 1 19861 ; iMathieu et all 
1990). This suggests that the relative velocity at infinity, and 
therefore the typical impact velocities of collisions, should be 
small compared to the significant orbital energy of compo- 
nent A in S1082 (and any very hard binaries that went into 
the encounter). From this, we take A m ~ 0. 

(iii) We are now ready to obtain quantitative constraints 
for the formation of S1082. We can use Equation [5] since 
M67 is an old OC with a small central velocity dispersion so 
that the assumptions used to derive this equation are valid. 
At this point, we must consider the details of our chosen 
formation scenario more carefully in order to choose a set 
of initial conditions that will satisfy Equation [9] as well as 
reproduce the observed parameters of S1082. In doing so, we 
find that it is not possible to form S1082 with a single 2+3 
encounter. This is because 5 stars with masses m < mro 
are insufficient to form both an inner binary with a total 
mass ~ 4.1 Mq and an outer tertiary companion with a 
mass 1.7 Mq. That is, the total mass of the inner binary is 
larger than three times the mass of the MSTO so that its 
formation must have involved four or more stars. 

A single 3+3 encounter is therefore the most probable for- 
mation scenario for S1082. Although there are a number of 
ways we can reproduce the observed parameters of S1082 
with a 3+3 encounter, including the component masses, we 
must use Equation [5] to identify the most probable of these 
scenarios. To do this, we can solve for the initial orbital en- 
ergies of all binaries and/or triples going into the encounter, 
which will in turn constrain their initial semi-major axes and 
therefore cross sections for collision. In applying Equation [9] 
we are only concerned with the initial and final orbital ener- 
gies of the inner binaries of the triples. This is because the 
orbital energies of any outer triple companions will be sig- 
nificantly outweighed by the orbital energies of their inner 



binaries. It follows that the contribution from the outer or- 
bital energies of triples to the total energy of the encounter 
will be negligible. 

From Equation [9] we find that the formation of S1082 
should have involved at least one hard binary with |e<| ~ 
10 41 J whose components did not merge (with each other) 
during the encounter in order to account for the significant 
orbital energy of component A. This need not be the case, 
however, if one or more stars were ejected from the system 
with an escape velocity of > 100 km/s. Of course, this would 
require an increase in the total number of stars involved 
in the interaction and therefore a single encounter scenario 
involving 7 or more stars, and hence one or more quadruple 
systems. Although very few observational constraints for the 
fraction of quadruple systems in clusters exist, this seems 
unlikely. 

The orbital energies of the inner binaries of the two triples 
initially going into the encounter should both be on the order 
of 10 41 J. Although we have found that only one hard binary 
with \ti\ ~ 10 41 J is required, triples are only dynamically 
stable if the ratio of their inner to outer orbital periods is 
large (roughly a factor of ten or more) fe.g. lMardlindl200ll ). 
We do not expect very wide binaries to survive for very long 
in dense cluster environments, which suggests that the inner 
binary of every triple is very hard. We expect from this that 
these inner binaries should, to within an order of magnitude, 
all have roughly the same orbital energy. 

The presence of outer triple companions is required in 
order for most very hard binaries to actually experience en- 
counters within the lifetime of a typical merger product. As- 
suming masses of 1 Mq for both components, an orbital en- 
ergy of 10 41 J corresponds to a period ~ 0.4 days, or a semi- 
major axis ~ 0.02 AU. Therefore, the cross section for colli- 
sion for a 2+2 encounter in which both binaries are very hard 
is smaller than the average cross section for a 2+2 encounter 
(found from the observed binary period distribution) by a 
factor of ~ 30. This suggests that the time required for an 
encounter to occur between two very hard binaries is consid- 
erably longer than the cluster lifetime. This is not the case 
if the hard binaries have triple companions, however, since 
the outer orbit significantly increases the collisional cross 
section and hence decreases the encounter time-scale. 

Energy conservation also suggests that if S1082 did form 
from a single 3+3 encounter, it is likely that the close inner 
binaries of the two triples collided directly. If two separate 
mergers then subsequently occurred during this interaction, 
this could have reproduced the observed orbital parameters 
of the close inner binary of S1082 (component A). The for- 
mation of the outer triple companion is more difficult to 
explain via a single 3+3 encounter since it also appears to 
be the product of a stellar merger. Nonetheless, if the outer 
companions of both triples undergoing the encounter end 
up orbiting the interacting (or merged) pair of close inner 
binaries at comparable distances, it is possible that their or- 
bits would overlap as a result of gravitational focussing. Al- 
though this seems unlikely, it could produce a blue straggler 
of the right mass and orbital period to account for compo- 
nent B. We will return to this possibility in Section [4] 

Even though the analytic estimates presented here are 
only approximate, they serve to highlight the low probabil- 
ity of a system such as S1082 having formed dynamically in 
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M67. Based on our results, the most likely formation sce- 
nario for S1082 is a single 3+3 encounter, although we ex- 
pect very few, if any, to have occurred in the lifetime of a 
typical merger product. This need not be the case, however, 
provided M67 had a higher central density in the recent past 
since this would increase the total encounter frequency. This 
will be discussed further in Section [4] 



3.2 The Period-Eccentricity Distribution of the 

Blue Straggler Binary Population in NGC 188 

iMathieu fe Gelled l|20Qgh found 21 blue stragglers in the old 
open cluster NGC 188. Of these, 16 are known to have a 
binary companion. Orbital solutions have been found for 15 
of the se known BS binaries. From this, IMathieu fe Gellerl 
(2009) showed that the BS binary population in NGC 188 
has a curious period-eccentricity distribution, with all but 
3 having periods of ~ 1000 days. Of these three, two have 
periods of < 10 days (binaries 5078 and 7782). Interest- 
ingly, one of these short-period BS binaries has a non-zero 
eccentricity. The normal MS binary population, on the other 
hand, shows no sign of a period gap for 10 < P < 1000 days 
l|Mathieu fe Ge llcr 2009). We can apply the procedure out- 
lined in Section 12.21 to better understand how we expect 
mergers formed during dynamical encounters to contribute 
to the BS binary population in NGC 188. Although the 
method described in Section ^. 21 treats one system at a time, 
we will apply our technique to the BS binary population of 
NGC 188 as a whole. 

(i) Before applying our technique, we must satisfy our- 
selves that a dynamical origin is possible for a large fraction 
of the observed BS population. Several examples of evidence 
in favour of a dynamical origin exist. For one thing, most 
BSs in NGC 188 have be en found to have binary companions 
|Mathieu fe Gelled2009f ). This should not be the case if most 
BSs are the products of the coalescence of isolated binary 
systems. On the other hand, a binary companion should be 
expected if the BS s were formed from m ass transfer. How- 
ever, the results of IChen fe Hani (|2008b) could suggest that 
most of the BSs in NGC 188 were not formed from mass 
transfer alone, providing indirect evidence that many BSs 
were formed from mergers. A binary companion should also 
be expected if the coalesced binary had a tertiary compan- 
ion to begin with. Although the Kozai mechanism can act to 
decrease the orbital separation of the inner binary of a triple, 
an additional perturbation is often required in order to fully 
induce the binary to merge (|Perets fe Fabrvcky||2009T ). All 
of this suggests that dynamical interactions likely played at 
least some role in the formation of a significant fraction of 
the BS binary population in NGC 188. 

The first step of our technique is to find the mos t com - 
monly occurring type(s) of encounter (sV lCeller et all (I2008T ) 
found a completeness-corrected multiple star fraction fi, + 
ft ~ 0.27 out to a period of ~ 4000 days. This represents a 
lower limit since it does not include binaries with P > 4000 
days. Usi ng the same low er limit for the ratio ft/ ft ~ 0.1 
found by iLathaml l|2005h . Figure [T] suggests that 1+2 en- 
counters should currently dominate in NGC 188. 

Second, we must constrain the minimum number of stars 
required to form the observed systems. From this, we can 
show that a merger occurring during all but a 1+1 encounter 



could produce a BS in a binary since this is the only type 
of encounter that involves less than 3 stars. By dividing the 
total mass of an observed system by mro, we can find an 
estimate for the minimum number of stars required to have 
been involved in its formation. Since all but one BS binary 
contain only a single merger product with a mass ^ tuto, 
realistic dynamical formation scenarios for these systems re- 
quire 3 or more stars. Binary 7782, on the other hand, is 
thought to contain two BSs so that if its formation involved 
two separate mergers it must have required 4 or more stars. 

A more quantitative comparison of the different encounter 
types is required since we do not yet know if enough en- 
counters occurred in the last tbs years to account for all 16 
observed BS binaries. To do this, we require estimates from 
the observations for a few additional cluster parameters in 
order to calculate and compare the relative encounter time- 
scales. NGC 188 ha s been found to have a central density of 
po = 10 2 ' 2 L© pc~ 3 dSandauist et al.ll2003l) an d a core radius 
of r c ~ 1.3 pc l|Bonatto. Bica fe Santosll2005l ). We found an 
average stellar mass for the cluster of < m >~ 0.9 M©. This 
was done by determining the cluster luminosity function us- 
ing only those stars kno wn to be cluster mem bers from the 
proper-motion study of IPlatais et alj (12003 ). W ith this, a 
theoretical isochrone taken from lPols et al.l 1 1998h was used 
to determine the cluster mass function and the average stel- 
lar mass was calculated. 

A sufficient number of suitable dynamical interactions 
should have occurred in the last tbs years for the forma- 
tion of every BS binary in NGC 188 to have been directly 
mediated by the cluster dynamics. Assuming once again 
that tbs ~ 1.5 Gyrs, the encounter time-scales derived in 
Appendix [A] suggest that a minimum of nine, eight, eight, 
three and one 1+2, 2+2, 1+3, 2+3 and 3+3 encounters, re- 
spectively, occurred within the lifetime of a typical merger 
product. It follows that at least ~ 29 dynamical encoun- 
ters should have occurred in NGC 188 in the last tbs years. 
Of these 29 encounters, 12 should have involved triples. In 
deriv ing these time-scal es, we have taken ft + ft ~ 0.27 
from iGeller at al.l (120091 ) and have adopted the same ratio 
ft/ fb ~ 0.1 as found for M67. Consequently, these represent 
lower limits for and ft , so that our derived encounter rates 
are also lower limits. We have also assumed an average outer 
semi-major axis for triples of at ~ 3 AU (corresponding to 
a period of ~ 1000 days for a binary composed of two 1 M© 
stars) . 

(ii) The next step is to apply our energy conservation 
prescription to the observed BS population. This will al- 
low us to constrain the orbital energies of typical binaries 
or triples expected to form BS binaries via dynamical in- 
teractions. First, we must refer to the literature in order to 
obtain estimates for every term in Equation [5] From the ob- 
served BS period-eccentricity distribution in NGC 188, we 
know that most BS binaries have periods of ~ 1000 days 
(we will call these long-period binaries), although there are 
a couple with very short periods of ~ 10 days (which we will 
call short-period binaries) . This provides an estimate for the 
final orbital energies e^k of BS binaries formed during dy- 
namical interactions. Specifically, we find \ekk I ~ 10 39 J and 
~ 10 40 J for the final absolute orbital energies of the 
long- and short-period BS binaries, respectively. 

Every BS in NGC 188 with a high cluster membership 
probability has both radial and proper motion velocities 
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that, to within their respective uncertainties, are consistent 
with the observed central velocity dispersion of op = 0.41 
km s " 1 (jPlatais et al.ll2003l : [Geller et al.ll200sl : iGeller at all 
2009). From this, we assume that the final translational ki- 
netic energies of any binaries or triples formed from dynam- 
ical interactions will be negligible, or Tjj ~ 0. As for M67, 
we assume A m ~ since we expect low impact velocities for 
collisions as a result of the low velocity dispersion in NGC 
188. 

(iii) We are now equipped with estimates from the liter- 
ature that will allow us to obtain quantitative constraints 
for specific encounter scenarios. In particular, we can use 
Equation [5] to constrain the initial orbital energies of all bi- 
naries and/or triples going into an encounter. We can also 
constrain the specific details of interactions for which the 
encounter outcome reproduces the observed parameters of 
the BS binary. We will consider energy conservation sepa- 
rately for two different classes of BS binaries, namely short- 
and long-period. 

The short-period BS binaries have large (absolute) orbital 
energies. Equation [9] suggests that this energy requires that 
at least one hard binary was involved in the encounter. Al- 
ternatively, an encounter which involved only softer bina- 
ries must have resulted in at least one star escaping with 
a significant velocity (> 100 km/s). However, in order to 
produce a merger product in a short-period binary, encoun- 
ters involving triples are the most favoured. Binary-binary 
(and especially 1+2) encounters which involve at least one 
hard binary will have smaller cross sections for collision than 
those encounters involving wide binaries. But energy conser- 
vation suggests that 2+2 encounters involving wide binaries 
are unlikely to produce a merger product in a close binary. 
Stabl e triples, on the other hand, contain a hard inner bi- 
nary (|Mardlinsll200lh which will naturally account for the 
large orbital energy of the resulting BS binary. Stable triples 
also have a large cross section for collision because of the 
wide orbit of the tertiary companion. The times between 
1+2, 2+2 and 1+3 encounters are all comparable, suggest- 
ing that most encounters involving hard binaries are 1+3 
encounters. 

The short-period BS binaries could have formed from a 
direct stellar collision that occurred within a dynamical en- 
counter of a hard binary and another single or (hard) binary 
star (we call this Mechanism I). If at least one of the objects 
going into the encounter was a triple, then four or more 
stars were involved in the interaction. Therefore, if binaries 
5078 and 7782 were formed from this mechanism, they could 
possess triple companions with sufficiently long periods that 
they would have thus far evaded detection. This is consis- 
tent with the requirements for both conservation of energy 
and angular momentum. Interestingly, the presence of an 
outer triple companion could also contribute to hardening 
these BS binaries via Kozai cycles operating i n conjunction 
with tidal friction (|Fabrvckv &: Tremainell2007h . In this case, 
Equation [l] shows that tides can contribute to making a bi- 
nary's orbital energy more negative by depositing internal 
energy into the component stars. 

Although 2+2 encounters should also contribute to the ob- 
served BS population, most of these should occur between 
a short-period binary and a long-period binary. This is be- 
cause most encounters resulting in mergers involve very hard 
binaries and binaries with long periods have large cross sec- 



tions for collision and therefore short encounter times. Equa- 
tion [9] indicates that the minimum period of a BS binary 
formed during a 2+2 encounter is usually determined by 
the orbital energy of the softest binary going into the en- 
counter. This is because the most likely merger scenario is 
one for which the hard binary is driven to coalesce by impart- 
ing energy and an gular momentum to other stars involved 
in the interaction |Fregeau et al .1 12004}). Assuming most of 
this energy is imparted to only one of the stars causing it to 
escape the system, we can take A m in Equation[9]to be very 
small and the orbital energy of the left-over BS binary will 
be comparable to the orbital energy of the initial wide binary 
going into the encounter. As more energy is imparted to the 
star left bound to the merger product over the course of the 
interaction, its final orbital separation effectively increases. 
This in turn contributes to an increase in the final orbital 
energy of the left-over BS binary. Finally, from Equation lA8l 
the period of the initial wide binary should be relatively long 
(roughly > 1000 days) since the contribution from the very 
hard binary to the total cross section for collision is negligi- 
ble. Otherwise, the time required for such a 2+2 encounter 
to occur could exceed tbs- 

Now let us consider the long-period BS binaries in NGC 
188. Based on Equation [5] we expect encounters involving 
3 or more stars and only one very hard binary to typically 
produce long-period BS binaries if the hard binary is driven 
to merge during the encounter. This is because the hard bi- 
nary merges so that its significant (negative) orbital energy 
can only be re-distributed to the other stars by giving them 
positive energy. Since the only other orbits involved in the 
interaction are wide, the left-over BS binary should also have 
a wide orbit. Alternatively, BSs formed during interactions 
involving more than one very hard binary should be left in 
a wide binary provided enough energy is extracted from the 
orbit of the binary that merges. In this case, a significant 
fraction of this energy must be imparted to the other stars 
in order to counter-balance the significant orbital energy of 
the other very hard binary. If not, the other short-period 
binaries are required to either merge or be ejected from the 
system. Otherwise Equation [9] suggests that the left-over BS 
binary should be very hard. However, it is important to re- 
call that we are neglecting other non-dynamical mechanisms 
for energy extraction. We will return to this last point in 
Section g] 

In order to help us obtain more quantitative estimates for 
the long-period BS binaries, consider two additional mecha- 
nisms for mergers during dynamical interactions that involve 
both short- and long-period orbits. First, a merger can oc- 
cur if a sufficient amount of orbital energy is extracted from 
a hard binary by other interacting stars (called Mechanism 
Ha). Alternatively, a merger can occur if the encounter pro- 
gresses in such a way that the eccentricity of a hard binary 
becomes sufficiently increased that the stellar radii overlap, 
causing the stars to collide and merge (called Mechanism 
lib). In the case of Mechanism lib, most of the orbital en- 
ergy of the close binary must end up in the form of internal 
and gravitational binding energy in the merger remnant af- 
ter the majority of its orbital angular momentum has been 
redistributed to the other stars (and tides have extracted or- 
bital energy). In the case of Mechanism Ha, however, most 
of the orbital energy of the close binary must be imparted to 
one or more of the other interacting stars in the form of bulk 
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kinetic motion. Consequently, one or more stars are likely to 
obtain a positive total energy and escape the system. This 
need not necessarily be the case for 2+3 and 3+3 encoun- 
ters provided the second hardest binary orbit involved in the 
interaction has a sufficiently negative energy. 

Regardless of the type of encounter, Equation [9] shows 
that the extraction of orbital energy from a hard binary 
in stimulating it to merge should increase the final orbital 
period of a BS binary. To illustrate this, we will consider 
1+3 encounters since the predictions of energy conservation 
are nearly identical for the other encounter types of interest. 
Moreoever, we have shown that encounters with triples are 
the most likely to involve very hard binaries. Equation [9] can 
be re-written for 1+3 encounters: 

612,4 = 612,3 ~ /l2 X £i2, (15) 

where we have assumed that stars 1 and 2 comprise the ini- 
tial hard inner binary of the triple, star 3 is the initial outer 
triple companion and star 4 is the interloping single star. 
Stars 1 and 2 are assumed to merge during the encounter 
by exchanging energy and angular momentum with stars 3 
and 4. We further assume that enough energy is imparted 
to star 3 that it escapes the system. We let /12 represent 
the fraction of energy extracted from the orbital energy of 
the hard inner binary of the triple in the form of bulk ki- 
netic motion by star 4. Since the remaining orbital energy of 
the close inner binary will end up in the form of internal and 
gravitational binding energy in the merger remnant, we have 
assumed (1 + /12) x £12 ~ T3 — A m in obtaining Equation ll5l 
from Equation [5] 

This formation mechanism could leave the BS as a single 
star if /12 exceeds a few percent. As shown in Figure [2l 
the period of a BS binary formed during a 1+3 encounter 
via Mechanism Ha (Pi.2,4) is only slightly smaller than the 
period of the outer orbit of the triple initially going into the 
encounter (Pi2,3). This assumes, however, that no energy is 
exchanged between the hard inner binary and star 4 (i.e. 
fi2 = 0). The predictions from energy conservation for this 
case are therefore identical for Mechanism lib. In general, as 
the amount of energy extracted from the hard inner binary 
of the triple by star 4 increases, so too does the rate at which 
Pi2,4 increases with increasing Pi.2,3. If the amount of energy 
extracted is > 5% of the orbital energy of the initial inner 
binary of the triple, Pi2,4 becomes a very steeply increasing 
function of Pi2,3- 

Interestingly, the two general qualitative merger scenar- 
ios described above (Mechanisms I and II) naturally create 
a bi-modal period distribution similar to the period gap ob- 
served for the BS binaries if we assume that 1+3 encounters 
produced these objects. To illustrate this, Figure [3] shows 
a histogram of periods for 15 BS binaries formed during 
1+3 encounters via these two generic merger scenarios. In 
order to obtain this plot, we have used the observed period- 
eccentricity distribution for the regula r MS-M S binary pop- 
ulation in NGC 188 from lGeller at all (2009) to obtain esti- 
mates for the orbital energies of any binaries and/or triples 
going into encounters. Specifically, in order to obtain pe- 
riods for the outer orbits of triples undergoing encounters, 
we randomly sampled the regular period distribution, in- 
cluding only those binaries with periods satisfying 400 days 
< P < 4000 days. We have shown that the initial outer 
orbits of triples going into 1+3 encounters provide a rough 
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Figure 2. Plot showing the typical periods of BS binaries formed 
from 1+3 encounters in which the hard inner binary of the triple 
merges. As described in the text, the period of the BS binary 
formed during the interaction is denoted P12 4 and corresponds 
to the y-axis, whereas the period of the outer orbit of the triple ini- 
tially going into the encounter is denoted P12 3 and corresponds to 
the x-axis. The lower straight line corresponds to the case where 
no energy was extracted from the inner binary of the triple by 
star 4 (i.e. /12 = 0). As the amount of energy extracted increases, 
however, so too will the rate at which Pi2,4 increases with in- 
creasing Pi2,3- Cases where /12 = 0.005, /12 = 0.01, /12 = 0.02, 
/12 = 0.03 and /12 = 0.04 are shown as lines of increasing slope. 



lower limit for the periods of BS binaries formed via Mech- 
anism II. Therefore, any BS binaries formed in this way 
could only have been identified as binaries by radial velocity 
surveys if the triple going into the encounter had a period 
< 4000 days (since this corresponds to the current cut-off 
for detection). All triples are taken to have a ratio of 30 
between their inner and outer orbital semi-major axes. This 
ratio has been chosen to be arbitrarily large enough that the 
triples should be dynamically stable, however we will return 
to this assumption in Section [4] 

We will adopt a rat io based on the observations of 
iMathieu fc Gelled ([2009) for the fraction of outcomes that 
result in each of these two possible merger scenarios. In par- 
ticular, if we assume that the three BS binaries with P < 150 
days were formed via Mechanism I whereas the other 12 
were formed via Mechanism II (either Ha or lib since en- 
ergy conservation predicts similar periods for the left-over 
BS binaries), this would suggest that Mechanism II is ~ 4 
times more likely to occur than Mechanism I during any 
given 1+3 encounter. We will return to this assumption in 
Section @] 

As a result of the requirement for conservation of angular 
momentum, we would expect to see a large spread in the 
distribution of eccentricities for BSs in long-period binaries 
formed from encounters involving triples. This is because 
conservation of momentum requires that the total momen- 
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Figure 3. Histogram of the period distribution (in days) expected 
for BS binaries formed during 15 1+3 encounters. The parameter 
space assumed for the encounters is described in the text. The 
stars show the observed BS binary pe riod distribution in NGC 188 
taken from Mathicu & Gcllcr (2009), where each star represents 
a single BS binary. 



turn going into the encounter must be equal to the total 
momentum left over after the interaction. However, the to- 
tal initial momentum depends not only on the initial orbital 
eccentricities of any binaries or triples going into the en- 
counter, but also the relative orientations and trajectories 
of the colliding objects. Since the relative orientations and 
trajectories are random, the final eccentricities of the BS bi- 
naries can take on a range of values. In other words, we can- 
not predict the final distribution of momenta for BS binaries 
formed from dynamical encounters. However, the observed 
BS binaries in NGC 188 are observed to have a range of 
eccentricities and this is not inconsistent with a dynamical 
origin involving triples. 

BSs in short-period binaries formed via Mechanism I can 
also end up with just about any eccentricity immediately 
after the encounter for the same reasons outlined above. 
However, tidal effects become increasingly significant with 
decreasing orbital separation so that the hardest binaries 
should typically circularize the fastest. Although theoretical 
estimates for the rat e of tidal circularization are uncertain 
l|Meibom et al.ll2005l ). the circularizati on cut-off period is es - 
timated to be ~ 15 days in NGC 188 l|Mathieu et al.ll2004) . 
From this, it is entirely possible that recently formed BS 
binaries with P ~ 10 days have not yet become fully circu- 
larized. 



4 SUMMARY & DISCUSSION 

In this paper, we have presented a generalized analytic 
prescription for energy conservation during stellar encoun- 



ters. Our method can be used to identify the most prob- 
able dynamical formation scenario for an observed binary 
or triple system containing one or more merger products. 
We have shown that, using the observed orbital parame- 
ters of the system, the allowed initial orbital semi-major 
axes of any binary or triple systems involved in its forma- 
tion can be constrained. The initial semi-major axes of the 
orbits in turn provide an estimate for the collisional cross 
section and therefore the time-scale for the encounter to 
occur in its host cluster. In order to apply our technique, 
repeated spectroscopic measurements of the binary or triple 
system containing the merger product (s) are needed in order 
to obtain its orbital solution and systemic velocity. How- 
ever, the time-scales provided in Appendix [X] can still be 
applied if only the fraction of binary and triple stars are 
know n, which can be determined eith er spectroscopically 
(e.g. Mathicu ct al. 1990; Latham 2005) or photometrically 
(e.g. ban et al.lll99g ). 

As we have shown, consideration of the requirement 
for energy conservation is ideal for identifying trends dur- 
ing stellar encounters, whereas numerical scattering experi- 
ments can require hundreds or even thousands of simulation 
runs before patterns will emerge. Some of these trends in- 
clude: 

• At least one short-period binary is usually required in 
a dynamical interaction to produce another binary having 
a similarly short-period (provided no stars are ejected with 
escape velocities > 100 km/s). This is because the orbital 
energy of a short-period binary is sufficiently negative that 
it tends to considerably outweigh the other energy terms in 
Equation \T\ for most of the encounters that typically occur 
in the cores of globul ar and especially op en clusters. This 
has been confirmed bv lHurlev et alj (|2005t ). 

• Previous studies have found that in order for triples to 
remain stable for many dynamical times, the ratio of their in- 
ner to outer orbital periods m ust be relatively large (roughly 
a factor of ten or more) (e.g. lMardlin3l200ll ). Based on our 
results, this has two important corollaries for stellar mergers 
in dense cluster environments hosting a significant popula- 
tion of triples: 

(i) The longest-lived triples will contain very hard inner 
binaries with a large |e|. This is important since stel- 
lar radii are in general more likely to overlap and hence 
mergers to occur durin g resonant interactions involving 
very hard binaries (e.g. iFregeau et al]|2004l : iHurlev et al.l 
2005). We therefore expect stellar mergers to be common 
during encounters involving stable triple systems. 

(ii) The longest-lived triples will contain wide outer or- 
bits, creating a large cross section for collision. This sug- 
gests that the time-scale required for a stable triple system 
to encounter another object is typically short relative to 
the cluster age in dense environments. A significant frac- 
tion of encounters involving very hard binaries, and hence 
resulting in stellar mergers, will therefore involve triples 
in old open clusters such as M67 and NGC 188. 



Ivan den Berg et all (|200ll ) and ISandquist et all (|2003l ) 
suggest that back-to-back binary-binary encounters, or even 
a single 3+3 encounter, could have formed S1082. We have 
improved upon these previous studies by estimating time- 
scales required for possible dynamical formation scenarios to 
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occur. Since we have argued in Section [3] that the formation 
of S1082 must have involved at least 6 stars, it follows that 
only a 3+3 interaction could have reproduced the observed 
configuration via a single encounter. However, the derived 
3+3 encounter time-scale is sufficiently long that we expect 
very few, if any, 3+3 encounters occurred within the life- 
time of a typical merger product. Moreover, we have argued 
that the times for multiple encounters to occur are longer 
than the cluster age. Although we cannot rule out a dynam- 
ical origin for S1082, our results suggest that it is unlikely 
(provided the derived encounter time-scales were not signif- 
icantly higher in the recent past, which we will return to 
below). From this, it follows that a dynamical link between 
the close binary and third star is unlikely to exist. 

On the other hand, we have so far ignored the cluster 
evolution, and assumed that the currently observed cluster 
parameters have not changed in the last few Gyrs. N-body 
models suggest that the central density in M67 could have 
been si gnificantly higher in the recent past. Specifically, Fig- 
ure 5 of lHurlev et al. (|2005T l indicates that the presently ob- 
served central density in M67 could have been higher within 
the past Gyr by a factor of > 2. If this was indeed the case, 
our previous estimates for each of the different encounter 
frequencies should increase by a factor of ~ 4, so that a sig- 
nificant number of dynamical encounters involving single, 
binary and triple stars should have occurred in M67 within 
the last tbs years. It follows that a dynamical origin for 
S1082 is not unlikely if the central density in M67 was re- 
cently larger than its presently observed value by a factor of 
> 2. This also increases the probability that a scenario in- 
volving multiple encounters created S1082, although we have 
shown that such a scenario is still likely to have involved one 
or more triples. 

Based on the preceding arguments, S1082 offers an ex- 
cellent example of how observations of individual multiple 
star systems containing BSs can be used to directly con- 
strain the dynamical history of their host cluster. If a defini- 
tive dynamical link between components A and B is estab- 
lished, this would suggest that the central density in M67 
was higher in the last 1-2 Gyrs. This is also required in order 
for the cluster dynamics to have played a role in the forma- 
tion of a significant fraction of the observed BS population 
in M67. Based on the current density, the encounter time- 
scales are sufficiently long that too few encounters should 
have occurred in the last tbs years for mergers during dy- 
namical interactions to be a significant contributor to BS 
formation. 

We have obtained quantitative constraints for two 
generic channels for mergers during encounters involving 
triples - one in which a direct stellar collision occurs within 
a dynamical interaction of the hard inner binary of a triple 
and another single or (hard) binary star (Mechanism I) and 
one in which the hard inner binary of a triple is driven to 
coalesce by imparting energy and/or angular momentum to 
other stars involved in the interaction (Mechanism II). Our 
results suggest that these two general merger mechanisms 
could contribute to a bi-modal period distribution for BS 
binaries similar to that observed in NGC 188. These dual 
mechanisms predict a gap in period, with those BS bina- 
ries formed via Mechanism I having periods of a few to 
~ 100 days and those formed via Mechanism II having pe- 
riods closer to ~ 1000 days. Some 2+2, 2+3 and even 1+3 



encounters could involve orbits with periods in this range, 
and Equation [9] confirms that the final period of a BS bi- 
nary formed via Mechanism II will typically be determined 
by that of the second hardest binary orbit. Therefore, we 
might still expect some BS binaries to have periods that 
fall in the gap (100 days < P < 1000 days). Our results do 
indeed predict one such BS binary, as shown in Figure [3] 

A number of assumptions went into obtaining Figure O 
many of which were chosen specifically to reproduce the ob- 
served BS binary period distribution. Regardless, our as- 
sumptions were chosen to reflect encounter scenarios that 
are the most likely to result in mergers. These should in- 
volve triples with very hard inner binaries since these are the 
most likely to merge during encounters. The triples should 
also have outer companions on very wide orbits since these 
have the largest cross sections for collision. From this, we 
have assumed a ratio of 30 between the inner and outer 
semi-major axes of all triples. This ensures that all triples 
are dynamically stable. This also leads us to assume a min- 
imum period of 400 days for the outer orbits of triples so 
that the corresponding minimum period of their inner or- 
bits is not too small. These assumptions serve to show that 
encounters involving triples could produce both long-period 
and short-period BS binaries as well as a period gap. 

Our results predict that the short-period peak in Fig- 
ure [3] is at a slightly longer period than the observations 
suggest. If we decrease the assumed ratio between the in- 
ner and outer orbital separations of triples, the short-period 
peak will move to even longer periods. However, if one or 
more stars were ejected with a high escape velocity or the 
dissipative effects of tides are considered (which are expected 
to be the most significant for encounters involving hard bina- 
ries) this would move the short-period peak to even shorter 
periods. In order to obtain the desired agreement with the 
observations at the short-period end of the BS period dis- 
tribution, energy must have somehow been dissipated or re- 
moved from the hard inner binaries of the triples during (or 
even after) the encounter, or the inner binaries must have 
initially been even harder than we have assumed in obtain- 
ing Figure [3] 

Our results could suggest that the hard BS binaries in 
NGC 188 (binaries 5078 and 7782) may have outer triple 
companions, perhaps with sufficiently long periods that they 
would have thus far evaded detection. This is consistent with 
the requirement for energy conservation since the orbital en- 
ergy of the outer orbit of the triple is negligible compared 
to that of its inner binary. If binaries 5078 and 7782 do 
have outer triple companions, it is also possible that Kozai 
oscillations combined with tid al friction contr ibuted to de- 
creasing their orbital periods l|Eggletonl l200rj ) . Finally, the 
BS binaries in NGC 188 are observed to have a wide range 
of eccentricities, which we have argued is not inconsistent 
with a dynamical origin involving triples. 

We have assumed that Mechanism II is more likely to 
occur than Mechanism I during any given 1+3 encounter. 
This is a reasonable assumption since numerical scattering 
experiments of 1+2 and 2+2 encounters have shown that 
the coalescence of a hard binary is much more likely to oc- 
cur during an encounter than a direct collision between a 
hard binary and a n interacting single or binary star (e.g. 
iFregeau et ai] [2004 ) . The ratio we have assumed for the fre- 
quencies with which Mechanisms I and II occur was specifi- 
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cally chosen in order to reproduce the observed numbers of 
short- and long-period BS binaries. The important point to 
take away is that the observed BS binary period-eccentricity 
distribution offers a potential constraint on the fraction of 
encounters that result in different merger scenarios. 

Based on our results, Mechanism 11 must occur ~ 4 
times more often than Mechanism I in order to reproduce 
the observed BS period distribution from 1+3 encounters 
(or, equivalently, 2+3 encounters involving a very wide bi- 
nary and 2+2 encounters between a short-period binary 
and a long-period binary). This can be tested by perform- 
ing numerical scattering experiments of encounters involving 
triples. Therefore, our results highlight the need for simula- 
tions of 1+3, 2+3 and 3+3 encounters to be performed in 
order to better understand their expected contributions to 
BS populations in open and globular clusters. Once a pre- 
ferred encounter scenario has been identified for an observed 
binary or triple containing one or more BSs, numerical scat- 
tering experiments can be used to further constrain the con- 
ditions under which that scenario will occur (or to show 
that it cannot occur). We have demonstrated that a combi- 
nation of observational and analytic constraints can be used 
to isolate the parameter space relevant to the dynamical for- 
mation of an observed multiple star system (or population 
of star systems) containing one or more merger products. 
This will drastically narrow the relevant parameter space 
for numerical scattering experiments. 

We have i mprov ed u pon the resu lts of 

IPerets fc Fabrvckvl (|2009h and iMathieu fc Gelled (120091 1 
since we have shown that dynamical encounters involving 
triples could not only be contributing to the long-period BS 
binaries in NGC 188, but they could also be an important 
formation mechanism for short-period BS binaries and 
triples containing BSs. We have not ruled out mass transfer 
or Koz ai-induced mergers in t r iples (primordial or other- 
wise) jMathieu fc Gelled 120091 ; IPerets fc Fabrvckvl 120091 ). 
or even various combinations of different mechanisms, as 
contributing formation channels to the BS binary popula- 
tion in NGC 188. For instance, a 1+3 exchange interaction 
could stimulate a merger indirectly if the resulting angle of 
inclination between the inner and outer orbits of the triple 
exceeds ~ 39°, ultimately allowing the triple to evolve via 
the Kozai mechanism so that the eccentricity of the inner 
binary increase s while its period remains roughly constant 
|EggletorJl2006l ). 

There is evidence to suggest that mass transfer via 
Roche lobe over flow could play a role in the formation of 
at least some BSs. It is difficult to account for the near zero 
eccentricities of some of the long-period BS binaries without 
at least one episode of mass transfer having occurred. This 
is because none of the normal MS-MS binaries with simi- 
lar p eriods have such small eccentricities (|Mathieu fc Gellerl 
120091 ). On the other hand, it may not be unreasonable to 
expect that some collision products left in binaries undergo 
mass transfer since they are expected to expand adiabati- 
cally post-collision, and will sooner or later evolve to ascend 
the giant branch. As a result of conservation of energy and 
angular momentum, the mass transfer process will usually 
act to increase the orbital per iods of these binaries provided 
it is conservative (jibed fl99ll ). Interestingly, the cut-off pe- 
riod f or Roche lobe overflow is ~ 1000 days for low-mass 
stars (|Eggletonll200l ). which is in rough agreement with the 



long-period peak in the observed period-eccentricity distri- 
bution of the BS binary population in NGC 188. Therefore, 
mass transfer could also be contributing to the period gap 
observed for the BS binaries. 

According to the results of iGeller at al.l (2009), the 
number of giant-MS binaries with P < 1000 days is com- 
parable to the number of BS binaries (A. Geller, private 
communication). It is unlikely that every giant-MS binary 
will form a BS from mass transfer, however, suggesting that 
at most a few of the long-period BS binaries in NGC 188 
were formed via this mechanism. Finally, if the outer com- 
panion of a triple system evolves to over-fill its Roche lobe it 
could transfer mass to both of the components of the close 
inner binary. This mechanism could therefore also produce 
two BSs in a close binary, although it predicts the presence 
of an orbiting triple companion. For these reasons, a better 
understanding of triple evolution, as well as binary evolution 
in binaries containing merger products, is needed. 

The dissipational effects of tides tend to convert stars' 
bulk translational kinetic energies into internal or thermal 
energy within the stars, leading to an increase in the total 
gravitational binding energy of t he stellar configuration (e.g. 
iMcMillan. Hut fc Makinol 1 19901 ). By increasing the terms 
Uii in Equation [T] the initial orbital energies of any binaries 
going into an encounter can increase accordingly in order to 
conserve energy. A higher orbital energy corresponds to a 
larger semi-major axis and hence cross section for collision. 
This suggests that the derived encounter time-scales can be 
taken as upper limits in the limit that tidal dissipation is 
negligible. We expect tides to be particularly effective dur- 
ing encounters for which the total energy is very negative as 
a result of one or more very hard binaries being involved. 

We have argued in Section 12.21 that the average stel- 
lar mass is expected to be comparable to (but slightly less 
than) the mass of the MSTO in old OCs and low-mass GCs. 
We have also argued that most encounters will involve stars 
having masses slightly larger than the average stellar mass. 
We might therefore expect to find that a high proportion of 
merger products have masses that exceed that of the MSTO 
in very dynamically-evolved clusters that have lost a large 
fraction of their low-mass stars. Consequently, a larger num- 
ber of merger products could appear sufficiently bright to 
end up in the BS region of the cluster CMD in these clusters 
than in their less dynamically -evolved counterparts . This is 
con sistent with the result s of Knigge. Leigh fc Sills! (|2009h 
and iLeigh. Sills fc Kniggd' (|2009l ) who found that the num- 
ber of BSs in the cores of GCs scales sub-linearly with the 
core mass. In particular, since the cluster relaxation time 
increases with increasing cluster mass, it is the lowest mass 
GCs that should have lost the largest fraction of their low- 
mass stars. Therefore, if a larger fraction of merger products 
do indeed end up more massive than the MSTO in these 
clusters, this could be a contributing factor to the observed 
sub-linear dependence on core mass. It is also interesting 
to note that, since BSs are among the most massive cluster 
members and many are thought to have a binary compan- 
ion, BSs should be preferentially retained in clusters as they 
evolve dynamically compared to low-mass MS stars. This 
could also contribute to the observed sub-linear dependence 
on core mass. 

The exchange or conversion of energies that occurs dur- 
ing an encounter takes place over a finite period of time, so 
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it is important to specify whether or not the system has fully 
relaxed post-encounter when disc ussing the r emain ing stel- 
lar configuration. For one thing, ISills et al l (|200lh showed 
that, although collision products may be in hydrodynamic 
equilibrium, they are not in thermal equilibrium upon for- 
mation and so contract on a thermal time-scale. Simulations 
also suggest that most merg er products should be rapid rota- 
tors (|Sills et al.ll2003 . 12003 ). However, at least in the case of 
blue stragglers, this is rarely supported by the observations. 
Some mechanism for angular momentum loss must therefore 
be operating either during or after the merger takes place in 
order to spin down the remnant. The time-scale considered 
must also be sufficiently short that subsequent dynamical 
interactions are unlikely to have occurred since these could 
affect the total energy and momentum of the system. 

With this last point in mind, N-body simulations con- 
sidering BS formation have shown that after they are 
formed, BSs are often exch anged into other multiple star 
systems l|Hurlev et al l l2005h . This suggests that for mul- 
tiple star systems containing more than one BS, the BSs 
could have first been formed separately or in parallel, and 
then exchanged into their presently observed configuration. 
For the case of S1082, this would require at least 3 sepa- 
rate dynamical interactions. Given that the derived times 
between encounters are relatively long and the fact that the 
most likely formation scenario is usually that for which the 
number of encounters is minimized, the current state of M67 
suggests that the probability of S1082 having formed from 
a scenario involving 3 encounters is low. Conversely, the 
derived encounter time-scales in NGC 188 are sufficiently 
short that many of the BS binaries could have experienced a 
subsequent dynamical interaction after their formation. BSs 
tend to be more massive than normal MS stars, contributing 
to an increase in their gravitationally-focussed cross section 
for collision. This suggests that the encounter time-scale for 
multiple star systems containing BSs is slightly shorter than 
for otherwise identical systems composed only of normal MS 
stars. This contributes to a slight increase in the probabil- 
ity that a BS will experience an exchange encounter after 
it is formed. Interestingly, it could also contribute to an 
increase in the probability that a close binary containing 
two BSs will form during an encounter between two dif- 
ferent multiple star syste ms each containing their own BSs 
|Mathieu &: Gellen ((2009) ; R. Mathieu, private communica- 
tion). This is because it is the heaviest stars that will experi- 
ence the strongest gravitational focussing and are therefore 
the most likely to experience a close encounter, end up in a 
closely bound configuration, or even merge. 

Most exchange interactions will involve wide binaries 
for which the cross section for collision is large. Since wide 
binaries are typically relatively soft and the hardest binary 
involved in the interaction will usually determine the orbital 
energy of the left-over BS binary, most exchange interactions 
will leave the periods of BS binaries relatively unaffected. 
This need not be the case, of course, provided one or more 
stars are ejected from the system with a very high escape 
velocity. With these last points in mind, we have assumed 
throughout our analysis that all binaries and triples are dy- 
namically hard. This is a reasonable assumption since the 
hard-soft boundary corresponds to a period of ~ 10 6 days 
in both M67 and NGC 188, for which the cross section for 
collision is sufficiently large that we do not expect such bi- 



naries to survive for very long. Nonetheless, considerations 
such as these must be properly taken into account when iso- 
lating a preferred formation scenario and predicting the final 
distribution of energies. 

We have presented an analytic technique to constrain 
the dynamical origins of multiple star systems containing 
one or more BSs. Our results suggest that, in old open clus- 
ters, most dynamical interactions resulting in mergers in- 
volve triple stars. If most triples are formed dynamically, 
this could suggest that many stellar mergers are the culmi- 
nation of a hierarchical build-up of dynamical interactions. 
Consequently, this mechanism for BS formation should be 
properly included in future N-body simulations of cluster 
evolution. A better understanding of the interplay between 
the cluster dynamics and the internal evolution of triple sys- 
tems is needed in order to better understand the expected 
period distribution of BS binaries formed from triples. Sim- 
ulations will therefore need to track both the formation and 
destruction of triples as well as their internal evolution via 
Kozai cycles, stellar and binary evolution, etc. On the ob- 
servational front, our results highlight the need for a more 
detailed knowledge of binary and especially triple popula- 
tions in clusters. 



APPENDIX A: COLLISIONAL CROSS 
SECTIONS AND TIME-SCALES 

The gravitationally-focused cross sections for 1+1, 1+2, 
2+2, 1+3, 2 +3 and 3+3 col lisions can be found using Equa- 
tion 6 from iLeonardl (|l989h . Neglecting the first term and 
assuming that binary and triple stars are on average twice 
and three times as massive as single stars, respectively, this 
gives for the various collisional cross sections: 



Cl+l 



01+2 



02+2 



01+3 



8-KGmR 



3irGmab 



8nGmah 



4nGma t 



0"2+3 



5irGm(ab + a t ) 



03+3 



12nGmat 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



Values for the pericenters assumed for the various types of 
encounters are shown in Table IA11 where R is the aver- 
age stellar radius, is the average binary semi-major axis 
and a t is the average semi-major axis of the outer orbits of 
triples. 

In general, the time between each of the different en- 
counter types can be found using Equation \W\ and the 
gravitationally-focused cross sections for collision given 
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Table Al. Pericenters Assumed for Each Encounter Type 



Encounter Type 


Pericenter 


1+1 


2R 


1+2 


a b /2 


1+3 


at/2 


2+2 


a b 


2+3 


(a b + a t )/2 


3+3 


at 



above. Following the derivation of iLeonardl (|l989h . we can 
write the encounter rate in the general form: 



^ X-{-y N X Tly(J X -l-yV X + y , 



(A7) 



where N x and n y are the number and number density, re- 
spectively, of single, binary or triple stars and v x + y is the 
relative velocity at infinity between objects x and y. For 
instance, the time between binary-binary encounters in the 
core of a cluster is given by: 



. 2+ ^1.3xl0V b - 2 (^) 3 (^) 
V r c / V no / 

/ Vrms \ / O.5M \ i \AU \ 

\5km/s J V< m > I V a b I 



(A8) 



] years. 



r c / \ no 

Vrms \ / O.5M \ t \AU \_ 

v 5km /s/\<m>/ \ a b 

Similarly, the times between 1+1, 1+2, 1+3, 2+3 and 3+3 
encounters are given by: 

. 1+1 ^i.ixio- ( i-/ b -/ t r 2 (^) 3 (^) 2 

V r c I V no / 



\ / 0.5Mq \ / 0.5-Rq \ 
V 5km/ s )\<m>)\< R>) 



years, 



(A9) 



ri+2 = 3.4 x 10 (1 - f b - ft) f b 



/ Vrms \ / 0.5M Q \ / 1AU \ 

\5km/ s ) V < m > / V a b ) 



V r c I \ no 
years, 



(A10) 



n+3 = 2.6 x 10 (1 - f b - ft) f t 

V r c 1 V n n / 

/ v rms \ / O.5M \ ( lA U \ 

J\<m>)\ at / 



V 5km/ s ) V < m > / V a t 

^^.oxioVrVr 1 ^)^^ 

JW_\ / 0.5M Q \ / \AU 
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and 
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